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\Q , Abstract 

, Suzaku has. for the first time, enabled the hard X-ray variability of the Seyfert 1 galaxy MGG-6-30-15 

' to be measured. The variability in the 14-45 keV band, which is dominated by a strong reflection hump, 

is quenched relative to that at a few keV. This directly demonstrates that the whole reflection spectrum is 
Q ■ much less variable than the power-law continuum. The broadband spectral variability can be decomposed 

I into two components - a highly variable power-law and constant reflection - as previously inferred from 

C/3 . other observations in the 2-10 keV band. The strong reflection and high iron abundance give rise to a 

' strong broad iron line, which requires the inner disc radius to be at about 2 gravitational radii. Our results 

are consistent with the predictions of the light bending model which invokes the very strong gravitational 
effects expected very close to a rapidly spinning black hole. 

Key words: galaxies: individual (MGG-6-30-15); galaxies: active; galaxies: Seyfert; X-rays: galaxies 



1. Introduction be about 2-3 times the Solar value (Fabian et al 2002), 

making the iron line particularly strong. 

The Seyfert 1 galaxy MCG-6-30-15 at z = 0.00775 has MCG-6-30-15 was observed several more times with 

been the centre of much interest since a broad iron line ASCA (Iwasawa et al 1996, 1999; Matsumoto et al 2003; 

was discovered in its X-ray spectrum with the ^S'Cvl satel- Shih et al 2002), by BeppoSAX (Guainazzi et al 1999), 

lite (Tanaka et al 1995). Iron line emission is part of the RXTE (Lee et al 1999; Vaughan & Edelson 2001) and 

reflection spectrum produced by the hard X-ray power- by XMM-Newton (Wilms et al 2001; Fabian et al 2002, 

law continuum in the source irradiating the accretion disc 2003; Vaughan & Fabian 2004). All of these observations 

(Guilbert & Rees 1988; Lightman & White 1988; Ross & confirmed the presence and general broad shape of the 

Fabian 1993) and is broadened and skewed by Doppler, iron line. Evidence that the black hole in MGG-6-30-15 

gravitational redshift effects, and light aberration and is rapidly spinning has been obtained from the extreme 

bending (Fabian et al 1989). The low energy extent of breadth of the line by Iwasawa et al (1996), Dabrowski 

the line can reveal the inner radius of the accretion disc et al (1997), Wilms et al (2001), Fabian et al (2002) and 

and thus the black hole spin (for reviews see Fabian et most recently Reynolds et al (2005), who determine a spin 

al 2000; Reynolds & Nowak 2003 and Fabian & Miniutti parameter a = 0.989. 

2006). The iron abundance in MCG-6-30-15 appears to The XMM-Newton work emphasizes that much of the 

^i^ 7. I , emission arise from smaller radii, between 2 — 6 (where 

e-mail: minmttiOast. cam. ac.uk „ ' s \ 

Tg = GM/c ). This may help to explain the otherwise puz- 
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zling spectral variations shown by the source. Shih et al 
(2002), Fabian et al (2002), Fabian & Vaughan (2003) and 
Taylor, Uttley & McHardy (2003) found that the spectral 
variability can be explained by a simple two component 
model consisting of a highly variable power-law contin- 
uum (of almost fixed spectral slope) and a much less vari- 
able reflection spectrum, which of course includes the iron 
line. The refiection (and iron line) strength do not follow 
the power-law continuum intensity, as expected in a sim- 
ple reflection picture. The small radius of much of the 
emission can however explain this behaviour when it is 
recalled that as well as strong gravitational red-shifting 
occurring in this region, there is also strong gravitational 
light bending. This can disconnect variations in the con- 
tinuum from those of the reflection and has led to the de- 
velopment of the light-bending model (Fabian & Vaughan 
2003; Miniutti et al 2003, 2004), which is a generalization 
of earlier work on the strong field regime (Martocchia & 
Matt 1996; Reynolds & Begclman 1997). 

If the source of the continuum emission (assumed to be 
an isotropic emitter) changes location close to the black 
hole, then, even if the continuum has a constant intrinsic 
luminosity, it appears to the outside observer to change in 
brightness. This is just due to gravity bending the light 
rays out of the line of sight by different amounts depend- 
ing upon the precise location of the source. Much of the 
radiation is bent down onto the disc, so the observed re- 
flection intensity changes little. The two-component vari- 
ability pattern and its light bending interpretation have 
recently found application in many other accreting black 
hole sources (Fabian et al 2004, 2005; Miniutti, Fabian & 
Miller 2004; Ponti et al 2006). 

One striking feature of the model reflection spectrum 
is a large refiection hump peaking at 20-40 kcV. This 
is where the disc albedo is highest; at lower energies 
the albedo is reduced by photoelectric absorption and 
at higher energies it is reduced (at a given energy) by 
Compton down-scattering of the photons (see e.g. George 
&; Fabian 1991). The presence of the Compton hump 
has been confirmed in the spectrum of MCG-6-30-15 by 
BeppoSAX and RXTE observations (Guainazzi et al 1999; 
Lee et al 2000; Fabian et al 2002). What those observa- 
tions have not done is to determine the variability of the 
Compton hump and show that it varies in the same way 
as the rest of the reflection spectrum^ . Here we carry out 
this important step with the Hard X-ray Detector (HXD) 
on Suzaku. 

We concentrate on the Suzaku data above 3 keV since 
there is a complex warm absorber in MCG-6-30-15 (Otani 
et al 1996; Lee et al 2001; Turner et al 2003,2004). The 
difference spectrum between low and high states of the 
source shows that absorption is minimal above 3 keV 
(Fabian et al 2002; Turner et al 2004) and High Energy 
Transmission Grating (HETG) spectra from Chandra 

^ The possible anticorrclation between the reflection continuum 
and iron line reported by Lee et al (2000) may in part be due 
to confusion caused by variation in the fitted power-law index 
(see discussion of a spurious reflection-photon index correlation 
found with RXTE data by Vaughan & Edelson 2001). 



show no absorption feature around 6.5 keV which would 
indicate absorption by species of intermediate ionization 
which could affect our fits (Young et al 2005). 

2. The Suzaku observations 

MCG-6-30-15 was observed four times by Suzaku 
(Mitsuda et al 2006), once between August 17-19 in 
2005 for about 45 ks and three times in January 2006 
with longer exposures. Here we focus on the three 2006 
observations performed between 9-14 (150 ks), 23-26 
(99 ks)and 27-30 (97 ks) January 2006. We use event 
files from revision 0.7 of the Suzaku pipeline. Version 
0.7 processing is an internal processing applied to the 
Suzaku data obtained during the Suzaku Working Group 
phase, for the purpose of establishing the detector calibra- 
tion as quickly as possible. Some processes that are not 
critical for most of the initial calibration and scientific 
studies, e.g., aspect correction, fine tuning of the event 
time tagging of the XIS data, are skipped in version 
processing, making the quality of the products limited 
in these directions compared with the official data sup- 
plied to guest observers. The XIS data were screened 
with XSELECT as standard (see e.g. Koyama et al 2006). 
The XIS products were extracted from circular regions 
of 4.3' radius centred on the source, while background 
products were extracted from two smaller circular regions 
offset from the source (and avoiding the chip corners with 
calibration sources) with a total area equal to that of the 
source region. The latest response and ancillary response 
files provided by the instrument teams were used. For 
the HXD/PIN (Takahashi et al 2006), instrumental back- 
ground spectra were extracted from time dependent mod- 
els provided by the HXD instrument team, based upon 
a database of non X-ray background observations made 
by the PIN diode to date. Since the background model- 
ing is the key issue for the hard X-ray measurement with 
the HXD, the HXD team has provided two independent 
background models, model A and B, which use different 
algorithms (Kokubun et al. 2006). Spectral analysis of the 
source spectrum using the two models was found to give 
statistically indistinguishable results for the three 2006 
January MCG-6-30-15 observations over the whole 14- 
45 keV band used here, and we use the background model 
A in our analysis. The response files appropriate for the 
XIS nominal position observation were chosen dated as of 
2006/08/07. 

We have extracted products for the three front- 
illuminated CCD XIS detectors (XISO, XIS2, and XIS3) 
and for the back-illuminated CCD (XISl). The XIS2 and 
XIS3 detectors are found to produce remarkably similar 
spectra in the whole band used here. The XISO spec- 
trum is slightly fiatter than that from the XIS2 and XIS3 
detectors so we proceed by co-adding just the XIS2 and 
XIS3 products in our analysis (see Yaqoob et al 2006 for 
mention of structures in the XISO and XISl spectra). In 
Fig. 1, we show the broad-band 0.5-12 keV background 
subtracted light curve from the XIS2 and XIS3 detectors 
during the three pointed observations in 2006. We also 
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Fig. 1. The light curve from the XIS2/XIS3 detectors in the 0.5—12 keV band. We also show the count rates chosen to select the 
High Flux (HF) and Low Flux (LF) states. They are chosen so that the HF and LF spectra have approximately the same number 
of counts (6.5 X lO''^ in the whole 0.5-12 keV band). 



show, as a reference, the count rate levels selected to de- 
fine the High Flux (HF) and Low Flux (LF) states which 
will be used to study the spectral variability of MCG- 
6-30-15. As normal for this source (thought to harbour 
a black hole with a mass of 3 x 10® il/0, McHardy ct 
al 2005) the light curve exhibits large amplitude and rel- 
atively short timescale variability with variations up to 
factors 2-3 in a few ks. 

3. The 3 12 keV XIS spectrum 

We start our analysis of the Suzaku data by considering 
the 3-12 keV time-averaged co-added spectrum from the 
XIS2 and XIS3 front-illuminated CCD detectors. (XIS re- 
sponse and ancillary files 20060213. rmf and 20060415. arf, 
with a 6 mm extract radius, were used.) The most im- 
portant feature in this energy band is the strong, skewed, 
relativistic Fe Ka line which characterizes the X-ray spec- 
trum of MCG-6-30-15 enabling us to explore the nature 
and geometry of the accretion flow close to the central 
black hole with much higher accuracy than in any other 
object so far. In Fig. 2 we show the ratio of the data to 
a simple power law model fitted in the 3-4 keV and 7.5- 
12 keV band and absorbed by the Galactic column density 
(4.08 X 10^" cm~^). The residuals clearly show the asym- 
metric and broad Fe Ka line profile (top panel). In the 
bottom panel, we superimpose the Fe Ka profile as ob- 
served with the XMM~Newton EPIC-pn camera in 2001. 
The agreement between the two instruments is remarkable 
and demonstrates the excellent level of the Suzaku XIS 
calibration even in the early stages of the mission. The 
2-10 keV flux is 4.1 x 10"" erg cm^^ g-i ^^le XMM~ 
Newton observation and 4.0 x 10"^^ erg cm"^ in the 
Suzaku one. Notice the absorption/emission structures in 
the blue wing of the relativistic line which are clearly de- 
tected by both the EPIC-pn and the XIS detectors with 
excellent agreement. 



3.1. The Fe K band and the relativistic line 

As a first attempt to fit the Fe Ka line profile we con- 
sider the simplest possible spectral model comprising a 
power law continuum absorbed by a column of neutral 
matter (fixed at the Galactic value Nh = 4.08 x 10-^" cm~^) 
and a set of Gaussian emission lines. We consider 
first three Gaussian emission lines: a narrow unresolved 
~6.4 keV emission line (the narrow component of the Fe 
line from distant matter), two narrow unresolved '^6.7 keV 
and ~6.97 keV absorption lines (the Fexxv and Fexxvi 
resonant absorption line already detected by Chandra, 
Young et al. 2005), and an additional Gaussian emission 
line with width free to vary (representing the broad rela- 
tivistic Fe hue). We obtain a fit with = 2550 for 2235 
degrees of freedom (dof). The narrow component of the Fe 
hue is at 6.43t°:°^ keV (Uncertainties are quoted through- 
out the paper at the 90 per cent confidence level) and has 
an equivalent width (EW) of only 30 ± 5 eV. Absorption 
lines are detected at 6.74 ±0.03 keV and 7.04 ±0.05 keV. 
The broad Fe hue is at 5.88 ± 0.05 keV, has a width of 
840 ± 40 eV and an EW of 305 ± 20 eV. 

However, the best-fitting model leaves clear residuals 
in the 3-12 keV band. In particular the broad Fe line 
is not properly modelled and the residuals suggest the 
introduction of a double Gaussian model with one broad 
Gaussian around 6.4 keV and an even broader one at lower 
energies to model the extended red wing. We obtain a 
significant improvement of the statistics with = 2398 
for 2232 (Table 1). The energies of the three unresolved 
lines are 6.41^g'Q2 keV for the narrow Fe emission line, 
and 7.04±0.05 keV, and 6.73±0.04 keV for the absorption 
lines. Their equivalent widths are 25±5 eV, — (12±8), and 
— (15 ±6) eV respectively. The upper limit on the EW of 
the narrow Fe Ka emission line (only 30 eV) indicates that 
reflection from distant matter plays a minor and marginal 
role in MCG-6-30-16, as already demonstrated by high- 
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Fig. 2. Top: The Suzaku XIS data divided by a power law 
model fitted in the 3—4 keV and 7.5-10 kcV band are shown 
in the most relevant 3—8 keV band. Bottom: The data from 
the XMM-Newton EPIC-pn detector (red) are superimposed 
on the XIS data (black). In both cases the model is a power 
law absorbed by the Galactic column density and fitted in the 
3-4 keV and 7.5-12 keV band. 

resolution spectroscopy with the HETG Chandra gratings 
(Lee et al 2002; Young ct al 2005). The two absorption line 
energies are slightly higher than the rest-frame energies of 
Fexxv and Fexxvi resonant absorption and are consistent 
with an origin in a common outflow with velocity of few 
thousand km s^^. In particular, they are both consistent 
in equivalent width and outflow velocity with the values 
measured with Chandra (2.0lg g x 10^ km s~^, Young et 
al 2005). The photoelectric absorption implied by such a 
thin highly ionized absorber is small and has a negligible 
effect on the fits described here. 

As for the broad Fe line, the parameters of the two 
broad Gaussian lines describing its profile are reported in 
Table 1 (first model "Double Gaussian"). They are both 
clearly resolved and their cumulative EW is 320 ± 45 eV, 
although this value should be taken with care since the 
continuum is a simple absorbed power law and does not 
include the reflection continuum which must be associated 
with the broad Fe line. It is interesting to note that the 
width of the Fe line core at ^ 6.45 keV is ct > 260 eV (i.e. 
FWIIM~ 31950 km s^^) suggesting that it originates in 



the outer disc at < 80 from the central black hole, while 
the red wing has a width > 760 eV indicating that the 
emitting matter is located within 6.5 rg from the centre, 
already implying, even by means of a very simple and 
phenomenological model, that the black hole in MCG-6- 
30-15 is most likely a spinning Kerr black hole in which 
the accretion disc extends down within the marginal stable 
orbit for a non-rotating Schwarzschild black hole (6 rg). 

3.2. A self-consistent reflection model 

The multiple Gaussian fit described above is merely a 
phenomenological description of the hard spectrum. The 
Gaussian emission lines and, most importantly, the broad 
Fe line are the clear signature of X-ray reflection and the 
above spectral model did not include any reflection con- 
tinuum. Here we build a much more self-consistent model 
in which the broad Fe line is computed together with the 
associated reflection continuum. We use a grid of models 
from Ross & Fabian (2005) to obtain the X-ray reflection 
spectrum. The X-ray reflection model has the spectral 
slope of the illuminating power law F, the ionization pa- 
rameter ^, the Fe abundance, and the normalization as 
free parameters. However, we forced the photon index to 
be the same as the power law continuum for consistency. 
Since the reflection model does not include Ni as an ele- 
ment, we also include a Ni Ka line with energy fixed at 
7.47 keV. 

In order to reproduce the relativistic broad Fe line pro- 
file, the reflection spectrum is convolved with a relativistic 
kernel derived from the Laor (1991) code. The relativis- 
tic blurring parameters are the emissivity indexes gin and 
gout (where the emissivity is e = r~'^'" within the innermost 
6 rg and £ = 7--9out outside), the inner disc radius rin, and 
the observer inclination i. The outer disc radius is fixed 
at its maximum allowed value of 400 rg. The choice of 
a broken power law emissivity profile is motivated by the 
previous long XMM-Newton observation of MCG-6-30-15 
(Vaughan & Fabian 2004). 

The model is applied to the co-added XIS2 and XIS3 
time-averaged spectrum in the 3-12 keV band and we 
keep the three unresolved Gaussian lines as above. The 
results of the spectral fitting are reported in Table 1 
( "Blurred Reflection" model) . We obtain a better fit with 
= 2360 for 2230 dof. The better fit with respect to 
the phenomenological Double Gaussian model described 
above is due to a better description of the overall rela- 
tivistic Fe line profile. In addition, the model accounts 
for some high-energy residuals above about 10 kcV which 
were present in the above best-fitting solution due to the 
unmodelled reflection component. 

We measure a relatively standard F = 2.18to;o6 '^O'^" 
tinuum slope which compares very well with previous re- 
sults with XMM-Newton. The reflection component con- 
tributes significantly to the hard spectrum of MCG-6-30- 
15 and we measure a reflection fraction of i? = 2.8 ± 0.9 
{R=l corresponds to the level of reflection expected from 
2tt sr), fully consistent with previous results from a si- 
multaneous XMM-Newton and BeppoSAX observation 
(Vaughan & Fabian 2004). The precise value of the re- 



Suzaku observations of MCG-6-30-15 



5 



Table 1. Results of spectral fits to the 3-12 keV X1S2 and XIS3 time-averaged co— added spectrum with the different models used 
to describe the relativistic Fe line of MCG-6-30-15. We present results for a phcnomcnological double— Gaussian fit and for a much 
more self— consistent relativistically blurred refiection model. For the reflection model, we measure an Fe abundance of 2.0'^^'^ times 
solar. We only report results for the relativistic Fe line here. A more detailed fit is presented in Table 2, where high-energy data 
from the HXD/PIN detector are also included. A subscript p indicates that the parameter reached its min/max allowed value. 

Double Gaussian 



Continuum Ka Red Wing Ka Blue Core /dof 

T E a EW E a EW 

1.96 ±0.02 5.38 ±0.10 840t™ 130 ±15 6.45+^;°^ 290 ±30 190 ±30 2398/2232 



Blurred Reflection 



Continuum Relativistic Blurring Reflector X^/dof 

r i Tin qin gout R ^ 

2.18;^!^^ 38 ±4 1.6to:3e5 4.6j:|^:^ 2.6 ±0.3 2.8 ±0.9 65 ±45 2360/2230 



flection fraction will be better constrained in a subsequent 
analysis when the high energy data from the HXD/PIN 
detector are considered as well. The reflector is only 
mildly ionized with ^ = 65 ± 45 erg cm and the in- 
clination is constrained to be 38° ±4°. 

As for the relevant relativistic blurring parameters, the 
XIS data are able to constrain the inner disc radius to 
be smaller than Tin < 2.2 r^. If, as it is customary, 
Tin is identified with the innermost stable circular orbit 
around a Kerr black hole, such a small 90 per cent up- 
per limit implies that the black hole spin in MCG-6-30-15 
is a > 0.917, i.e. the black hole is an almost maximally 
spinning Kerr black hole. The inner and outer emissiv- 
ity profiles are not consistent with each other within the 
errors and gin c:i 4.8 is steeper than gout — 2.6. However, 
if the two indices are forced to be one and the same, we 
obtain gin = gout = 3.1^q3 with only a slight worsening of 
the fit statistic (Ax^ = 6). 

4. Adding the HXD/PIN data; the 3 45 keV 
spectrum 

In the above analysis the relativistic Fe line is associated 
with its own reflection continuum and a self-consistent 
model to the 3-12 keV spectrum is found. However, one of 
the main characteristics of the Suzaku mission is the pres- 
ence of the HXD PIN hard X-ray detector providing high 
quality data above 12 keV (Fig. 3). The spectral analysis 
of the XIS data revealed the presence of a strong reflection 
component associated with the broad relativistic Fe line 
which would produce a strong Compton hump around 20- 
30 kcV. The HXD data above 12 keV are thus crucial to 
investigate the refiection component further and will en- 
able us to infer the reflection parameters and to measure 
the direct and reflection continua with high accuracy. 

In Fig. 4, we show the XIS and HXD/PIN time- 




so 
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Fig. 3. PIN data displayed as total countratc (blue, upper), 
backgrounds A and B (red and green, respectively, middle) 
and source (black lower) spectra. 

averaged spectrum in the 0.5-45 keV band. The Figure 
shows the data to model ratio for a power law model fitted 
above 3 keV and ignoring the bands where reflection is ex- 
pected to dominate. The ratio clearly shows the presence 
of both the relativistic Fe line in the XIS data and also the 
presence of a large Compton hump around 20 keV visu- 
ally confirming that the hard spectrum of MCG-6-30-15 is 
dominated by a strong reflection component. Absorption 
is present affecting the soft band below 3 kcV. We avoid 
the complication of modelling the complex warm absorber 
and we concentrate, in this first paper on the Suzaku ob- 
servation of MCG-6-30-15, on the data above 3 keV. A 
detailed model of the warm absorber will be presented 
elsewhere where the data from the back-illuminated (more 
sensitive at soft energies) XIS detector will be also con- 
sidered. 

We consider our best-fitting spectral model to the 3- 
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Fig. 4. The XIS2/XIS3 and PIN data in the 0.5-45 keV band 
are compared with a simple power law fitted in the 3-45 keV 
band ignoring the Fe K and Compton hump energy bands 
(4-7.5 keV and 14-30 keV). The ratio plot shows the main 
spectral components which are the power law continuum, ab- 
sorption below about 3 keV, and a strong smeared reflec- 
tion component comprising the relativistic broad Fe line and 
Compton hump. Data have been re-binned for visual clarity. 

12 keV XIS data and extend it to the PIN data between 
14 keV and 45 keV. We include a cross-normalization 
constant between the XIS and PIN data and find that 
it is constrained to be 1.10 ±0.05 (where the XIS con- 
stant is set to 1). This factor is consistent with the 
results obtained from Crab observations (Kokubun et 
al. 2006). Since the available HXD/PIN background 
model A does not include the contribution from the X- 
ray cosmic background, a spectral model of the form 
2.06 X lQ-^{E/im fcel^)-i-29exp(-i;/41.13 keV) was in- 
cluded in all spectral fits in the HXD/PIN band (the 
model being valid up to about 70 kcV). The model is 
based on the HEAO-Al spectrum, re-scaled to account 
for the HXD field of view. 

The best-fitting XIS model provides an acceptable de- 
scription of the 3-45 keV data and we obtain a good qual- 
ity fit of = 2495 for 2312 dof. The model parameters 
are very much consistent with those obtained from the 
3-12 keV parameters. In fact, extending the model up 
to 45 keV provides a reasonable description of the hard 
data, reproducing most of the large Compton hump seen 
in Fig. 4, even before fitting. The results for the Time- 
Averaged (TA) spectrum in the 3-45 keV band are re- 
ported in Table 2 together with similar fits to the High 
Flux (HF) and Low Flux (LF) spectra which will be dis- 
cussed below. The refiection component contributes sig- 
nificantly to the total flux and we measure a large re- 
flection fraction of 3.8 ± 0.7. As a reference, and since 
it will be useful in a spectral variability analysis below, 
we report here that its contribution represents the 33 ± 6 
per cent of the total flux in the 3-12 keV band and the 
51 ± 10 per cent in the 14-45 keV band. The relativistic 
blurring parameters are consistent within the errors with 
those already reported in Table 1 and the same is true 
for the continuum power law slope and for the ionization 



Fig. 5. The best-fit model components to the time-averaged 
spectrum of MCG-6-30-15. 

parameter. 

As a consistency check, we have also determined the 
level of the reflection component by using the spectrum 
above 7.5 keV, so avoiding the iron line. For this the 
pexrav model in XSPEC was used, with results F = 2.2, 
reflection fraction i? = 2 — 4 and iron abundance poorly 
constrained (Apo = 0.5 — 2 times solar). We note that 
such reflection predicts a strong iron line with an equiva- 
lent width of at least 300 eV (taking conservatively an iron 
abundance of unity and reflection fraction of 2) , using the 
predictions of George & Fabian (1991). The limit on any 
narrow component to the iron emission feature is at least 
ten times smaller than this (this work; Lee et al 2002; 
Young et al 2005), thereby ruling out any non-smeared 
origin for the reflection component. The reflection must 
be due to matter well within the optical broad line re- 
gion. We also stress that if the evidence for relativistic 
reflection is taken into account and the pexrav compo- 
nent is blurred accordingly, the shape of the reflection 
continuum is broader and redshifted requiring a higher 
Fe abundance (^pc = 1.5 — 3.5 times solar) and reflection 
fraction (i? = 3 — 5) to fit the data, while the continuum 
photon index is basically unaffected. The blurred and Fe 
overabundant reflector described by pexrav is thus fully 
consistent with our best-fitting results with a much more 
complex spectral model including the relativistic Fe line 
(see Table 2). 

5. Flux and spectral variability 

The HXD/PIN detector has no imaging capabilities 
and thus there is no simple way to obtain a background 
subtracted light curve by selecting appropriate regions. 
However, background subtracted light curves can be ob- 
tained by extracting spectra in each time interval, sub- 
tracting the background from the available models in the 
same time interval, and reading out the background sub- 
tracted count rate. In background model A, the back- 
ground of the HXD is calculated from the data base by 
using the count rate of the upper discriminator, which has 
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Fig. 6. The background subtracted light curve of 
MCG-6-30-15 during the three 2006 pointed observa- 
tions (first in the top panel, second and third in the bottom 
panel). The dotted black light curve is the XIS2 and XIS3 
light curve rescaled to the mean HXD/PIN count rate in 
each orbit. In red, we show the HXD/PIN light curve. The 
time bin is about 45 ks. Time elapses from the start of the 
first observation. 

a strong correlation with the flux of the non- X-ray back- 
ground. The procedure must also include a model for 
the cosmic X-ray background (not included in the avail- 
able background models) and we use the same model as 
for spectral fitting by simulating the cosmic X-ray back- 
ground HXD/PIN spectrum and adding it to the instru- 
mental background. We have applied the above procedure 
to the HXD/PIN data by selecting, after some experimen- 
tation, a timescale of 45 ks enabling us to obtain about 
5000 background subtracted counts per time interval in 
the 14-45 keV PIN band on average. In Fig. 6 we show 
the XIS2 and XIS3 and the HXD/PIN light curves for 
the three 2006 observations in the 0.5-12 keV and 14- 
45 keV respectively. To ease the comparison, the XIS 
light curves have been re-scaled so that they have the 
same mean count rate as the PIN ones in each segment 
(i.e. observation). It is visually clear that the harder band 
(solid red) has a much lower variability amplitude than the 
softer one (dotted black). In the following we explore a 
set of different techniques with the aim of exploring the 



Fig. 7. Top: The XIS2 and XIS3 count rate in the hard 
3-12 keV band is plotted against that in the 1—2 keV band 
from the 45 ks light curve. The relationship is linear and 
reveals a clear hard off— set. Bottom: Flux— flux plot for 
the 14-45 keV HXD/PIN data against the 1-2 keV reference 
band. Due to larger spread in the PIN band (see text for 
details), we have re-binned the original flux-flux plot accord- 
ing to the 1-2 keV flux so that each data point shown here 
comprises three original data points. The binned flux- flux 
plot exhibits a linear relationship, although some scatter is 
present. 

flux and spectral variability of the source. 
5. 1 . Flux-flux plots 

As a first model-independent way to characterize the 
flux/spectral variability of MCG-6-30-15 we analysed the 
relationship of the fluxes (count rates) in different energy 
bands following the technique of flux-flux plots introduced 
by Churazov, Gilfanov & Revnivtsev (2001) in a study 
of Cygnus X-1 spectral variability, and Taylor, Uttlcy & 
McHardy (2003) as a model-independent tool to disentan- 
gle the main drivers of the spectral variability in AGN. 

In Fig. 7 we show the 3-12 keV XIS count rate from the 
45 ks light curve against the count rate in the 1-2 keV 
band which is chosen as reference (see also Vaughan & 
Fabian 2004). The relationship between the count rates 
in the two bands is remarkably linear and leaves a clear 
hard offset on the y-axis, as already noticed by Taylor, 
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Uttley & McHardy (2003) and Vaughan & Fabian (2004) 
in previous RXTE and XMM-Newton observations. The 
hncarity of the relationship indicates that the flux varia- 
tions arc strongly dominated by changes in the normal- 
ization of a spectral component with constant spectral 
shape. On the other hand, the hard offset suggests the 
presence of a component that varies little and that con- 
tributes more to the hard than to the reference 1-2 keV 
band. We measure a hard offset of a = 0.48 ±0.02 cts/s 
which can be used to infer the contribution of the weakly 
variable component at mean flux level: since the mean 
count rate in the 3-12 kcV band is 1.79 ±0.01 cts/s, the 
hard offset a represents 26.8l}'2 P^r cent of the 3-12 keV 
flux at mean flux level. This compares very well with the 
33 ± 6 per cent contribution of the reflection component 
to the 3-12 keV band, as obtained from direct spectral 
fitting. 

We have applied the same technique to the HXD/PIN 
light curve and plotted the 14-45 keV count rate against 
the 1-2 keV reference one. The relationship is much more 
noisy than when data from the XIS alone were used. To 
reduce the scatter and decrease the PIN error bars, we 
have binned the original flux-flux plot according to the 
soft reference count rate so that each new point comprises 
three original data points. In the lower panel of Fig. 7 
we show the binned flux-flux plot obtained by plotting 
the 14-45 keV HXD/PIN count rate against the reference 
1-2 keV one. 

A fit with a linear relationship is perfectly acceptable 
(x^ = 5.8 for 5 dof), again indicating the presence of a clear 
hard offset in the PIN band. In this case, the contribution 
of the weakly variable (reflection) component to the 14- 
45 kcV band at mean flux level is estimated to be as high 
as 55 ± 17 per cent. Again, this is in good agreement with 
the contribution of the reflection component in the 14- 
45 keV band obtained from direct spectral fitting (51 ± 10 
per cent). 

The flux-flux plot analysis detailed above indicates that 
the spectral variability of MCG-6-30-15 can be decom- 
posed into two main components: a highly variable com- 
ponent which varies in normalization only but not in spec- 
tral shape, and a weakly variable one which has a much 
harder spectral shape. The excellent agreement between 
the contribution of the weakly variable component in the 
different energy bands as inferred from the flux-flux plot 
analysis and that of the reflection component as obtained 
by direct spectral fltting strongly suggest that the weakly 
variable component is the rclativistically smeared X-ray 
reflection from the accretion disc. Obviously, the highly 
variable component can be identifled with the other com- 
ponent required by spectral fltting, i.e. the power law 
continuum. The flux-flux plots suggest that as the source 
varies, the power law maintains the same spectral slope 
and only varies in normalization. 

5.2. The Flux RMS spectrum 

A further approach to studying the spectral variabil- 
ity of a source is the RMS spectrum, which, for a chosen 
timescale, is the fractional variability seen at each energy 




a 5 10 20 

Energy (keV) 



Fig. 8. The RMS spectrum of MCG-6-30-15 on a 45 ks 
timescale is shown in the 2-45 keV band. The solid red line 
represents the theoretical RMS that would be obtained from 
the best-fitting spectral model of the time— averaged spectrum 
if all the variability could be explained with a two-component 
model in which the power law has a constant photon index 
and varies in normalization only, while the disc reflection com- 
ponent is perfectly constant. 

expressed clS &- flux normalized root mean square vari- 
ability function. The RMS spectrum allows us to quan- 
tify the fractional variability as a function of energy in a 
model-independent manner and is therefore an important 
piece of information, complementary to the flux-flux plot 
analysis described above. Such RMS spectra have previ- 
ously been determined for MCG-6-30-15 (Matsumoto et 
al 2003; Fabian et al 2002; Ponti et al 2004). However, as 
far as data below 10 keV are available, the RMS spectrum 
has proved to be ambiguous to model. In particular, the 
typical trend that is observed in AGN is that the fractional 
variability decreases with energy, a behaviour that can be 
explained either in terms of a power law softening at high 
flux levels (i.e. a pivoting power law) or in the framework 
of the two-component model discussed above. The new 
ingredient here is to include the HXD/PIN data above 
10 keV with the goal of confirming/rejecting the interpre- 
tation that comes out from the fiux-flux plots analysis. 

Our result for a 45 ks timescale is shown in Fig. 8. The 
general trend is that the RMS decreases with energy with 
a marked drop at ~ 6.4 keV reassuringly conflrming that 
the Fe line is less variable than the continuum, as already 
noticed in previous works (e.g. Vaughan & Fabian 2004; 
Ponti et al 2004) . The novelty of the Suzaku observation 
is represented by the HXD/PIN data. They confirm the 
trend of lower fractional variability at higher energies but 
they lie above and not on the extrapolation of the 2-6 keV 
trend. 

We have theoretically modelled the RMS spectrum by 
considering our best-fitting spectral model to the time- 
averaged spectrum and reproducing the 45 ks data by al- 
lowing only the power law normalization to vary, i.e. en- 
forcing a perfect two-component model. The result of this 
exercise is shown as a solid red line in Fig. 8 and shows 
that the RMS spectrum is fully consistent with a two- 
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Fig. 9. Ratio plot for the HF (black) and LF (red) spectra 
with a power law model fitted ignoring the 4-7.5 kcV and 
14-25 keV energy bands. The ratio shows that the Compton 
hump contributes by about 80 per cent in the LF state and by 
about 40-50 per cent in the HF one, while its contribution in 
the time— averaged spectrum is about 60 per cent (see Fig. 4). 
This shows that the reflection component is less variable than 
the power law and thus, the reflection fraction is expected to 
increase as the overall flux decreases. 

component model in which the only variable component 
is a constant-slope power law. Some variability of the re- 
flection component cannot be excluded, but its variability 
amplitude is certainly not comparable to the variability 
amplitude of the power law. Indeed, small mismatches 
between the model and the data could be due to some 
intrinsic variability of the reflection component which de- 
serves more detailed study to be performed in future work. 

5.3. The High and Low Flux states 

Having analysed the spectral variability of MCG- 
6-30-15 with two model-independent (and calibration- 
independent) techniques we now consider a more direct 
approach through spectral fitting. We extracted High 
Flux (HF) and Low Flux (LF) from both the XIS and the 
PIN spectra according to the selection criterion shown in 
Fig. 1. Wc first consider a simple power law fit in the 
3-45 kcV band ignoring, as done before, the Fc K and 
Compton hump bands. The result for the HF and LF 
spectra is shown in Fig. 9 as a ratio plot. There is a clear 
evidence, both in the XIS and in the HXD/PIN, that the 
reflection component contributes more to the LF than to 
the HF spectrum, which again indicates an almost con- 
stant reflection. 

We have then applied the same spectral model as for 
the 3-45 keV time-averaged (TA) spectrum to the HF 
and LF spectra and our results are reported in Table 2. 
Almost all parameters are consistent with being the same 
in the HF and LF spectra. In particular, the power law 
slope is consistent with being constant and equal to ~ 2.25 
in the TA, HF, and LF spectra. The same is true for the 
Fe abundance, disc inclination, inner disc radius, outer- 
most emissivity index (jout and reflector ionization state 
(although there is some indication for a lower ionization 
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Fig. 10. The spectra, best-fltting model, and residuals for 
the HF (upper panel) and LF (lower panel) states. 

at low flux levels). The innermost emissivity index q^^ ap- 
pears to be steeper in the LF than in the HF spectrum, 
but the two indexes are still consistent with each other 
within the 90 per cent errors. A steeper emissivity profile 
in the LF state would produce a more extended red wing 
as expected in the framework of the light bending model 
where low flux states are characterized by a more centrally 
concentrated illumination of the inner disc (Miniutti & 
Fabian 2004). 

On the other hand, the reflection fraction definitely 
appears to vary and it is higher in the LF spectrum 
(i? = 4.8±0.8) than in the HF one (i? = 2.5±0.6) as already 
strongly suggested by the ratio plot in Fig. 9. This result 
is fully consistent with the indication inferred from the 
fiux-fiux plot and rms spectrum analysis discussed above: 
if the refiection component is only weakly variable, its 
contribution increases in the LF states because the direct 
power law decreases while the reflection stays approxi- 
mately constant. We point out that this is exactly the be- 
haviour predicted by the light bending model (Miniutti & 
Fabian 2004). The spectra, best-fitting models and resid- 
uals for the HF and LF states are shown in Fig. 10. The 
best-fit spectral components are as in Fig. 5 ( which refers 
to the best-fit parameters to the time-averaged spectrum) 
with the only major difference in the HF/LF state being 
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Fig. 11. The difference spectrum (HF minus LF) fitted with 
s power law absorbed by the Galactic column. 

due to a higher/lower power law normalization. 
5.4- The difference spectrum 

We also produced a difference spectrum obtained by 
subtracting the LF spectrum from the HF one for the 
XIS and HXD/PIN. By definition, the difference spec- 
trum only shows the variable components of the spectrum 
(modified by absorption), while any non variable compo- 
nent is subtracted away. In Fig. 11 we show the residuals 
to a simple power law fit in the 3-45 keV band where 
Galactic absorption is included. The power law fit is very 
good and the power law slope is consistent with that found 
in our spectral analysis (F = 2.2 ± 0.1) but clear positive 
residuals are left in a broad hump around 6 keV at the 
~ 2(7 level. This may indicate either that the red wing of 
the relativistic Fe line has varied in fiux or that the shape 
of the line is slightly different in the HF and LF spectra. 

If the residuals are due to flux variability of the reflec- 
tion component, we should also see some positive residuals 
at the Compton hump energy (20-30 keV) which are not 
seen in the HXD/PIN difference spectrum. Some reflec- 
tion variability is in fact allowed by the results of our spec- 
tral fits of the HF and LF spectra (see Table 2) although 
the refiection fiux is consistent with being the same at the 
two flux levels. If the refiection flux varied but the Fe line 
keeps the same shape, the main residual in the difference 
spectrum should be at the energy of the blue peak of the 
line, around 6.4 keV, while the residuals in Fig. 11 indicate 
that only the red wing of the line has varied. 

The remaining possibility is that the Fe line has a differ- 
ent shape (rather than a different flux) in the HF and LF 
spectra. This is already hinted at by the inner emissivity 
index (gin, sec Table 2) which is marginally steeper in the 
LF state (although consistent within the errors). If the 
emissivity is really different at the two flux levels, resid- 
uals are produced in the red wing of the line and could 
potentially explain those seen in the difference spectrum 
(and at the same time the lack of residual at Compton 
hump energies). The result is in the sense expected from 
the light-bending model (see Fig. 4 of Miniutti & Fabian 



2004). 

6. Discussion 

The Suzaku observation of MCG~6-30-15 has, for the 
first time, enabled the hard X-ray variability of the source 
to be determined. A wide range of techniques including 
direct spectral fitting, fiux-fiux plots and rms spectrum 
all show that the spectral variability on a 45 ks timescale 
can be decomposed into just two broad components above 
3 keV, namely a variable power-law continuum and a 
harder constant component. This last component has the 
broad iron line and refiection hump expected from X-ray 
reflection in the very innermost regions of an accretion 
disc around a rapidly spinning black hole. 

Just the strength and shape of the reflection component 
predicts a high equivalent width iron line, much stronger 
than the weak narrow emission component seen in the 
spectrum. The only way for such a strong iron line to be 
'hidden' in the spectrum is for it to be smeared into the 
broad emission structure that we observe. 

The behaviour of the power-law and reflection compo- 
nents is consistent with the light bending model. In which 
case much of the rapid (power-law component) variabil- 
ity seen in the source is due to strong gravitational light 
bending as the site of the emission region changes loca- 
tion within a few gravitational radii of the black hole. We 
assume that the emission region is associated with some 
magnetic structure produced by the strong differential mo- 
tions in the accretion disc. Note that to produce observed 
variability the emission region does not need to physically 
move, just that the main emission site changes position. 

The broad iron line is strong in MCG-6-30-15 because 
both the iron abundance is high, at around 2 times the 
solar value, and the large degree of light bending around 
its Kerr black hole has emphasised the reflection spectrum 
relative to the power-law continuum. 
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Table 2. Results of fits to the XIS2/XIS3 and PIN Time-Averaged (TA), High Flux (HF) and Low Flux (LF) spectra in 
the 3-45 keV band. The power law normalization (A'^pl) is given in units of 10~^ ph cm~^ s~^. As for the reflection, 
we prefer to give the reflection fraction which is a more interesting measure of its strength. The power law and reflection 
fluxes are reported in the 3-45 keV band in units of 10~^^ erg cm~^ s~^ for completeness. The reflector ionization param- 
eter (^) is in erg cm and the inner disc radius (fin) and inclination (i) are in units of rg = GMjc? and degrees respec- 
tively. The energies of the lines are in keV, their EW in eV, and their width is fixed to 10 eV. The narrow Fe Ka nor- 
malization is in units of 10^® ph cm^-^ s~^. A subscript p indicates that the parameter reached its min/max allowed value. 



Parameter 


TA spectrum 


HF spectrum 


LF spectrum 


r 

A^PL 


2.26 ±0.04 
1.7±0.1 


2.25 ±0.05 
2.7±0.1 


2.25 ±0.05 
1.2±0.1 


R 

Tin 

i 

Q'out 


3.8 ±0.7 

68 ±31 

^ - ' -0.465p 

38 ±3 
4.4t°;i 
2.5 ±0.3 


2.5 ±0.6 
70 ±37 

37±4 

o n+0.8 

O.D_g 7 

2.8 ±0.4 


4.8 ±0.8 

48 ±35 
1 q+i-5 

1 7-1-0.3 
^- ' -0.465p 
41 ±5 

4.5 ±0.6 

2.6 ±0.3 


£^FcKq 
EWpcKa 


6 41+0 03 

35 ±5 
8.5±1.3 


6.40 ±0.03 
19 ±4 
8.3±1.2 


6.41 ±0.03 
44 ±6 
8.5±1.3 


EW,bs, 


6.73 ±0.04 
-(22 ±8) 
7.07 ±0.04 

-(13±8) 


6.72 ±0.05 

-(201?^) 
7.08 ±0.05 
-(15±10) 


6.75 ±0.04 
-(45 ±30) 
7.03 ±0.04 
-(28 ±18) 


pPL 

-^3-45 koV 
pREF 
-^3-45 kcV 


4.0 ±0.2 

4.1 ±0.8 


6.4±0.2 
4.0 ±0.9 


2.9 ±0.3 
3.8 ±0.6 


xVdof 


2495/2312 


1910/1848 


1900/1925 
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